Quantum Phase Transition between (Luttinger) Liquid and Gas of Cold Molecules 
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We consider cold polar molecules confined in a helical optical lattice similar to those used in holo- 
graphic microfabrication. An external electric field polarizes molecules along the axis of the helix. 
The large-distance inter-molecular dipolar interaction is attractive but the short-scale interaction is 
repulsive due to geometric constraints and thus prevents collapse. The interaction strength depends 
on the electric field. We show that a zero-temperature second-order liquid-gas transition occurs at 
, a critical field. It can be observed under experimentally accessible conditions. 

o ■ 
o ■ 

' At zero temperature most substances exist in the form of solids. The only exception is Helium which undergoes 
quantum melting at a critical pressure. On the other hand, zero-temperature liquid-gas transitions have never been 
' observed. Indeed, at absolute zero any system must be in its ground state and condensed phases have lower energy 
^**^ , than gases. Driving a quantum liquid into a zero-temperature gas state would be possible, if one could control 
inter-atomic forces. At sufficiently weak inter-atomic interactions, a many-particle bound state (i.e., liquid or solid) 
would cease to exist and a gas would form instead. While this cannot be accomplished with conventual materials, 
recent progress in the field of cold dilute gases opens a possibility to tailor a wide range of Hamiltonians with tunable 
parameters. In this Letter we show that a quantum liquid-gas transition can be observed in a cold gas of polar 
molecules confined in an optical lattice. 

Experiments with cold gases have already allowed the observation of Bose-Einstein condensation, BCS superfluidity 
and Mott localization^. It was proposed that cold gases can serve as realizations of other analogies of electronic matter 
such as superconductors with p-wave pairing 2 and quantum Hall states 3 . Besides, several new states of matter with 
different broken symmetries and/or soft modes were predicted in cold atom systems. We address a rather different 
situation. Liquids and gases have the same symmetry and their only difference is the density: A gas fills all available 
volume while the density of a liquid is determined by inter- molecular interactions. 
. While cold gases do not represent true ground states of alkali metals, they are highly stable due to their low 
density which guarantees a low probability of multi-particle recombination processes^. A well-established way to 
control interaction in cold gases utilizes Feshbach resonances and allows changing both the strength and sign of the 
short-range inter-atomic potential which can be modeled as a delta-function in space. If it is repulsive the system 
minimizes its energy by occupying all available volume, i.e., is a gas. Attractive interactions have different effects on 
£NJ ■ fermions and bosons. A Fermi gas undergoes Cooper pair formation while a Bose gas collapses into a regime with 
strong many-body recombination. In our case a different inter-atomic potential in needed: A liquid state can be 
formed, if the interaction is attractive at large scales but the short-range force must be repulsive to prevent collapse. 

We demonstrate that such potential can be built from the dipole interaction of polar molecules^. The sign of the 
dipole interaction is certainly independent of the distance and depends only on the direction of the dipole moments. 
• ■ It was shown recently^ that the interaction sign can be made distance-dependent by driving molecules with microwave 
fields. However, such an approach generates many-body forces^. It would be interesting to investigate if a liquid phase 
is possible in such a system but in this paper we focus on a simpler situation with only two-body interactions. It can 
be achieved by confining polar molecules in a helical optical lattice, i.e., a potential well of a helical shape. Hexagonal 
arrays of such helices^ are among numerous periodic and aperiodic structures used in holographic microfabrication 
experiments. We will see below that such structures open new possibilities for tailoring cold-atom Hamiltonians which 
k> cannot be obtained with usual optical lattices in the form of sinusoidal waves. 

We focus on the simplest problem of this type with the lattice in the form of a single helix (Fig. fa). It can 
be produced with an approach similar to Refs. as discussed in the Appendix. An external electric field polarizes 
molecules along the axis of the helix. In zero field, molecules have zero angular momentum and hence zero average 
dipole moment while in strong fields the average dipole moment p can reach values of the order of Debyes. At large 
inter- molecular distances the dipole interaction V = p 2 (l — 3 cos 2 6) /r 3 is attractive (6 = 0). At short distances it 
becomes repulsive, if the angle 7 between the tangent to the helix and its axis exceeds the magic angle cos _1 (l/\/3). 
The distance dependence of the interaction may include multiple maxima and minima for 7 w tt/2. In this paper we 
focus on smaller 7 so that the interaction dependence on the distance s along the helix has a simpler shape shown in 
Fig. lb. The potential resembles the Lennard-Jones interaction used in models of thermal liquid-gas transitions and 
we show that a transition between a Luttinger liquid and Fermi or Tonks-Girardeau gas occurs at a critical value of 
the dipole moment p. 

The paper is organized as follows: we first formulate the model. Next, we find its phase diagram with a variational 
method. We support the variational calculation by a proof that for a class of inter- molecular potentials there is a phase 
transition between a monoatomic gas and a liquid. Then we address the conditions for the experimental observation 
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FIG. 1: a) Helical lattice, b) Potential V(s) in units of p 2 /(nR) 3 as a function of distance s in units of nR. 7 = arctan2. 
of the transition. 

We consider N 3> 1 particles of mass m confined in a helical potential well of radius R and pitch d ~ i?. We 
assume that the wave function is confined in a region of the width ~ w <C i? around the helix (Fig. la). The adiabatic 
approximation applies and at low temperatures T <C h 2 /(mw 2 ) the problem reduces to a one-dimensional model with 
the Hamiltonian 

where Sj is the length of the helix between particle number i and a reference point on the helix, V(s) the dipole 
interaction. Repulsive short-range interaction keeps particles apart. Hence, in ID the particle statistics is unimportant. 
In what follows we focus on non-identical particles with s„+i > s n for all n. The ground state energy is independent 
of statistics and the wave functions of identical bosons and fcrmions can be obtained from the case of non-identical 
particles by symmetrization or antisymmetrization. 

As the first step, we use the self-consistent harmonic approximation (SCHA), well-known in the theory of disordered 
classical systems 8,9 . A better approximation will be considered in the end of the paper. Note that a method, similar 
to SCHA, has been used in a related problem^. We interpret the Euclidean action S as a classical Hamiltonian and 
h as the temperature. We consider a trial distribution function of the form P ~ exp(— Sq/K), where the trial action 
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contains the variational parameters h and K; h is related to the average length of the TV-particle chain and K 
describes the fluctuations of the particle positions Sk- We need to minimize the variational free energy E{K,h) = 
[{S — So)o — h\nZ ]/T 7 where (. . .) denotes the average with respect to the trial distribution function P, the trial 
partition function Zq = J Tldsk exp(— Sq/K), and the time r — > 00 is the size of the integration domain. The self- 
consistent harmonic approximation provides analytic results for potentials V(s) in the form of the sums of exponentials. 
Thus, we approximate the potential represented in Fig. lb by the Morse potential^, V(s) = A[exp(— 2a[s — ho}) — 
2 exp(— a[s — ho])], which has a similar shape. Since the interaction rapidly decreases at large distances, we keep only 
the interaction between neighboring particles (fc and k + 1). 

The variational calculation includes three steps: 1) we change the variables Sk — > Sfc + kh so that h drops 
out of Eq. @; 2) calculate E(K,h); and 3) find the minimum of E(K,h). Step 2) requires the calcula- 
tion of the correlation function^ (s ra (w)s^(aj)) = hG~^ n (uj)/2. The matrix G nm (ui) differers from a tridiago- 
nal Toeplitz matrbsi^ only by the values of two matrix elements. Its inverse can be found analytically using 
the same method as for tridiagonal Toeplitz matrices 1 -. In the large N limit, step 2) yields E(K,h)/[Nr] = 

, where x = exp(— a[h — ho]). After the minimization 

with respect to x, the energy per particle E(K,h)/[Nr] reduces to e(K) = (h/n)yj (K/m) — Aexp[—2a 2 h/TT^/Km}. 
The energy e(K) has a local minimum at K = 0, where e(0) = 0. At large A another minimum at K = K m i n is 
possible. The phase transition into the state corresponding to that minimum occurs when e(K m i n ) — 0. Solving the 
latter equation together with de(K m \ n ) / dK = we find a phase transition at A c — 2ea 2 h 2 /[mn 2 ]. 

Thus, at small A we get K = and h — 00, i.e. the system is a gas. K and h are finite at large A. This means a 
finite volume at zero external pressure, i.e. a condensed state. In ID this cannot be a solid and the calculation of the 
correlation function ((s n+ k(t) — s n (t) — kh) 2 ) w ril " fc at k 3> 1 shows that the particles form a Luttinger liquid^. 
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Our problem is connected with the physics of atoms confined in Carbon nanotubes. If one ignores the periodic 
potential created by Carbon atoms then a model, similar to ours, emerges (certainly, one cannot tune the interaction 
between the atoms and obtain a phase transition in a nanotube). A variational studylS. predicted that an increasing 
interaction drives a monoatomic gas into a diatomic gas phase before a liquid state can be reached. This contradicts 
our findings. Below we sketch a proof that a liquid state has lower energy than a di- or multi-atomic gas and hence 
a multi-atomic gas cannot be the ground state. Qualitatively this reflects the fact that in a liquid every particle has 
bonds with two neighbors and only one bond in a diatomic molecule. Hence, one expects lower e nerg y per particle 
in a liquid. The prediction of a diatomic gas is thus an artifact of the variational method (and Ref. 10| admits such a 
possibility) . 

Our proof is based on a variational upper bound for the ground state energy. We will also use that bound to 
improve SCHA. We will focus on models in which only neighboring particles (fc and fc + 1) interact. A diatomic gas 
was predicted in Ref. [13 for such a model and such a model is relevant for us since dipole forces rapidly decrease 
at large scales. We will discuss elsewhere a derivation of the statements, proven below, for systems in which all 
pairs of particles interact via attractive potentials with hard-core repulsion. We first demonstrate that the energy of 
a triatomic gas is always lower than the energy of a diatomic gas, provided that the interaction potential and the 
ground state wave functions are well-behaved. Then we use a similar argument to show that the energy of the liquid 
phase is lower than the energies of all possible multi-atomic gases. In all cases we assume that the potential energy 
is zero at an infinite interparticle separation Sfc+i — Sk = +°o and impose the hard core condition V^(0) = oo. 

Consider a system of N particles in the interval — oo < s^. < oo. The center of mass is at rest in the ground 
state. Hence, the wave function can be represented as tPn(Ai, . . . , Ajv-i), where Afe = s^+i — s^. The ground state 
configuration can be viewed as a set of molecules. A molecule is a bound cluster of n particles fc, (fe + 1), . . . , (fc+n— 1) 
with finite^ interparticle distances A&, . . . , Ak+ n -2- The molecules are separated by infinite^ intervals A p . Different 
molecules do not interact. Hence, the ground state energy is the sum of the energies of separate molecules: One can 
easily see that if the ground state includes k molecules with N%, .. . , Nk particles then the energy equals the sum of 
the ground state energies of the Hamiltonians (TTJ) with N = N%, . . . , Nk- 

We now show that the energy of a diatomic gas can always be decreased, if particles rearrange into trimers. A 
diatomic gas may form, if two particles have a bound state with the energy £2 < 0. The bound state wave function 
-02 (Ai) is a normalized eigenfunction of the Hamiltonian H%2, given by Eq. U) with N = 2. Since the Hamiltonian 
H12 = —^-js 1 + ^(^i) i s re& l' we can assume that ^2 is reali 7 . In the diatomic gas, the energy per particle is 
£2/2. We now demonstrate that there is a three-particle state with the energy per particle e' < £2/2. The three- 
particle Hamiltonian i/123 = — Y^i=i ^~§^ + ^(Ai) + V"(A2). Let us find the average energy E of the trial wave 
function ip3( s ii s 2> S3) — V-'2(Ai)V'2(A2). The structure of the trial wave functions prompts the change of variables 
Si, S2, S3 — > Ai, A 2 , S3. The Jacobian of this transformation is one. Hence, 



The integral of the last term in the square brackets reduces to the square of the integral of a full derivative, 
(/ dA[dtp2l dA]/2) 2 — 0. The first four terms in the brackets can be represented as the sum of two two-particle 
Hamiltonians Hn + Hiz- Thus, E — 2e2 and the energy per particle E/3 < £2/2. This shows that the energy is lower 
in the triatomic gas than in the diatomic gas which thus cannot be the ground state. 

A similar argument proves a general statement: Consider a system of N particles in the interval — 00 < s^ < 00. 
Only nearest neighbors interact. Assume that a bound state exists for n < N particles. Then in the ground state 
the system consists of no more than n molecules and exactly one of those molecules contains more than one particle. 
Indeed, the ground state energy of the infinite system equals the sum of the energies of its molecules with zero center- 
of-mass velocities. Imagine that the ground state includes two multi-atomic molecules with m and (fc + 1) atoms. Let 
tp m (Ai, . . . , A m _i) and ifjk+i(Ai, . . . , A*,) be the ground states of the Hamiltonian fl} with N — m and N = fc + 1 
and the energies of these states be e m and ek+i- The sum of these two energies contributes to the ground state 
energy e g . We now use the argument of the previous paragraph to show that the energy decreases, if we substitute 
those two molecules with a monoatomic molecule (whose energy is 0) and a (m + fc)— atomic molecule. Indeed, the 
same calculation as in Eq. ((3|) with the variational wave function ifj m +k = V'm(Ai, . . . , A m -i)il)k+i{A m , . . . , A m +fc_i) 
shows that the ground state energy t m +k of the Hamiltonian H m+ k (HJ with N = m + fc cannot exceed e m + e^+i- 
The wave function ip m +k is not an eigenfunction of H m +k and hence e m +k < e m + £fe+i- Thus, there is a state whose 
energy e g — [e m + £fc+i) + (0 + £&+,„) is below the ground state energy e g . The contradiction means that no more 
than one multi-atomic molecule exists in the ground state. Hence, if there were more than n molecules in the ground 
state then at least n of them would be monoatomic. In such situation the energy decreases, if we form an additional 
n-atomic molecule from n free particles. This proves that there are no more than n molecules and no more than one 
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FIG. 2: Laser beam configuration. 



of them is multi-atomic. If N 3> n then exactly one molecule contains at least N — n + 1 particles which can be 
described as a liquid. 

The above discussion shows that only two possibilities exist for the ground state in the large- N limit: a monoatomic 
gas, if no bound states exist at all, and a liquid. Both possibilities take place at different interaction strengths A. 
At A — the system must be a gasi£. At large A, SCHA yields a negative upper bound for the energy. Hence, the 
monoatomic gas with its zero energy cannot be the ground state at large A and a liquid-gas transition must occur at 
an intermediate A. SCHA is exact at small and large A. Indeed, it correctly predicts zero energy in the gas phase at 
small A. At large A the fluctuations of A& are small. Hence, it is legitimate to expand the potential energy V(A) up 
to the second order which means that the action is quadratic and SCHA is quantitatively valid. However, SCHA is 
insufficient near the phase transition. This is clear from the comparison of the variational estimate for the transition 
point A c and the exact threshold Ad for the formation of diatomic molecules in the Morse potential^. Contrary to 
the above proof, the SCHA result for A c exceeds Ad — h 2 a 2 /(Am). Thus, a different method is needed near the 
transition. We try the variational ansatz of the form tpN = ^k=i ip(Ak), where the whole function ip is a variational 
parameter. From the calculation, completely analogous to Eq. ([3]), we find the variational energy E = (N — l)i?2, 
where E 2 is the average energy of a two-particle system in the state ip(si — s^) in the Morse potential. The lowest E 2 
corresponds to ^>(A) which is the ground state in the Morse potentiate. Hence, E 2 = — A[l — ah/ %/ Am A] 2 at A > Ad 
and E = at A < Ad- This improves an estimate for the transition point, A c ^ ncw = Ad- Obviously, E is lower than 
the SCHA variational energy e = in the interval A c > A > A c _ ncw . The improved variational method leads to an 
unexpected prediction concerning the order of the liquid-gas transition. The size of a diatomic molecule diverges^ as 
s ~ \/\[~E at A — > Ad- Hence, the size I of the iV-atomic bound state ipN diverges according to the same law. This 
means that the density of the liquid p = N/l ~ \[E — > at A — > A c . In other words, the variational method predicts 
a second order liquid-gas transition. Second order transitions between Luttinger liquid states with different densities 
are known 1 ^ but in contrast to Refs. ,19 the symmetry does not change at our transition. According to the Landau 
theory such transitions must be first-order. On the other hand, examples— are known of second order transitions 
in low-dimensional systems for which the Landau theory predicts the first order. It would be interesting to find a 
rigorous description of the ID liquid-gas critical point. 

We see that the phase transition occurs when the characteristic kinetic and potential energies are of the same order 
of magnitude, (ha) 2 /2m ~ V(\/a). For polar molecules in a helical lattice, the characteristic spatial scale 1/a ~ ttR. 
Hence, near the transition the dipole energy p 2 /(irR) 3 must be of the order of the recoil energy h 2 /m(irR) 2 . For a 
realistic optical lattice with ttR ~ 1/im and the molecular mass of the order of 100 we find p ~ ID at the transition. 
Such dipole moments are within reach. A difference between a liquid and gas can be detected in a variant of the 
Einstein's-boxes experiment. A laser beam, orthogonal to the helix, creates a potential barrier in the center of the 
helical lattice. A gas occupies all lattice and an approximately equal number of particles will remain on both sides 
of the barrier. The volume of a liquid is much smaller than the system size far from the transition. Hence, all atoms 
will be found on one side of the barrier. 

In conclusion, we have shown that a gas of polar molecules in a helical optical lattice can be driven by an electric 
field into a Luttinger liquid state via a continuous phase transition. The gas fills all available volume while the volume 
of the liquid is determined by the interaction strength. We thank G. P. Crawford, J. Dalibard, E. Dcmler, A. Kitaev 
and C. Salomon for useful discussions. We acknowledge the support by NSF under Grant No. DMR-0544116 and the 
hospitality of LPTENS (D. E. F.). 

Appendix. A helical lattice can be obtained as shown in Fig. 2. A circularly polarized wave with 
E c = E c (l,i,0) and the wave vector fc(0,0, ±1) interferes with n laser beams with wave vectors k m = 
fc(cos a cos[27rm/n], cosa sin[27rm/n], sina) and electric fields e m = E\(— sin[27rm/n], cos[2-7rm/n], 0). At n = 6 
this configuration produces a periodic array of identical helices 7 . We focus on the mathematically simplest case 
of large n. The optical potential perceived by an atom in the L z — state is proportional to the intensity of light 4 
| E | 2 = 2E 2 + [n E\ J i (kp cos a)] 2 + 2nE c E\J \(kp cos a) cos[fcz(±l — sina) + </>], where 0, p, z are polar coordinates. De- 
pending on the sign of detuning the atoms will be trapped near the intensity minimum or maximum. Both correspond 
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to helices fcz(±l — sin a) + 4> — nn; p — const. 
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